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Abstract 

Nowadays it has become feasible to perform angle-resolved photoemission spectroscopy (ARPES) 
measurements of transition-metal oxides with three-dimensional perovskite structures owing to 
the availability of high-quality single crystals of bulk and epitaxial thin films. In this article, 
we review recent experimental results and interpretation of ARPES data using empirical tight- 
binding band-structure calculations. Results are presented for SrVOs (SVO) bulk single crystals, 
and Lai„a,Sr2.Fe03 (LSFO) and La1_2.Sr2.MnO3 (LSMO) thin films. In the case of SVO, from 
comparison of the experimental results with calculated surface electronic structure, we concluded 
that the obtained band dispersions reflect the bulk electronic structure. The experimental band 
structures of LSFO and LSMO were analyzed assuming the G-type antiferromagnetic state and 
the ferromagnetic state, respectively. We also demonstrated that the intrinsic uncertainty of the 
electron momentum perpendicular to the crystal surface is important for the interpretation of the 
ARPES results of three-dimensional materials. 
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INTRODUCTION 



Perovskite-type 3d transition-metal (TM) oxides have been attracting much interest in 
these decades because of their intriguing physical properties, such as metal-insulator transi- 
tion (MIT), colossal magnetoresistance (CMR), and ordering of spin, charge, and orbitals jlf. 
Early photoemission studies combined with cluster-model analyses have revealed that their 
global electronic structure is characterized by the TM 3d-3d on-site Coulomb interaction en- 
ergy denoted by U, the charge-transfer energy from the occupied O 2p orbitals to the empty 
TM 3c? orbitals denoted by A, and the p-d hybridization strength denoted by (pda) Q]- If 
A > U, the band gap is of the d-d type, and the compound is called a Mott-Hubbard-type 
insulator. For A < U, the band gap is of the p-d type: the O 2p band is located between the 
upper and lower Hubbard bands, and the compound is called a charge-transfer-type insula- 
tor. In going from lighter to heavier 3d TM elements, the d level is lowered and therefore 
A is decreased, whereas U is increased. Thus Ti and V oxides become Mott-Hubbard-type 
compounds and Mn, Fe, Co, Ni, and Cu oxides charge-transfer-type compounds. 

Having established the global electronic structure picture of the transition-metal oxides 
with their chemical trend, the next step is to understand the momentum dependence of 
the electronic states, that is, the "band structure" of the hybridized p-d states influenced 
by the Coulomb interaction U. Angle-resolved photoemission spectroscopy (ARPES) is a 
unique and most powerful experimental technique by which one can directly determine the 
band structure of a material. However, ARPES studies have been largely limited to low- 
dimensional materials and there have been few studies on three-dimensional TM oxides with 
cubic perovskite structures because many of them do not have a cleavage plane. Recently, 
Yoshida et al. succeeded in cleaving single crystals of SrVOa (SVO) along the cubic (100) sur- 
face and performed ARPES measurements pj]. Another way of performing ARPES studies 



of such materials is the use of single-crystal thin films. Recently, high-quality perovskite- 
type oxide single- crystal thin films grown by the pulsed laser deposition (PLD) method have 
become available and setups have been developed for their in-situ photoemission mea- 

surement jsl Q]. By using well-defined surfaces of epitaxial thin films, we recently studied 
the band structures of Lai-a'SrajFeOa (LSFO) and La1_a.Sri.MnO3 (LSMO), and it has been 
demonstrated that in-situ ARPES measurements on such TM oxide films are one of the 
best methods to investigate the band structure of TM oxides with three-dimensional crystal 
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structures 



In this paper, we report on analyses of ARPES results of bulk gugle crystals of SVO Q 
and on m-sttu ARPES results of single-crystal LSFO (x = 0.4) j9| and LSMO (x = 0.4) 
0] thin films grown on SrTiOs (001) substrates. In the case of SVO, we have investigated 
the effect of surface on the ARPES spectra by comparing the experimental results with 
calculated surface electronic structures. The experimental band structures of LSFO and 
LSMO were interpreted using an empirical tight-binding band-structure calculation by as- 
suming the G-type antiferromagnetic state for LSFO and the ferromagnetic state for LSMO. 
The mass enhancement compared with band-structure calculation is discussed in the case 
of SVO and LSMO. We also discuss the problem of intrinsic uncertainty in determining the 
momentum of electrons in the solid perpendicular to the crystal surface jlo| . 

THEORETICAL ANALYSES 

Tight-binding band-structure calculation 

In the case of SVO, we adopted the simplest tight-binding (TB) model. In this simplest 
model, the t2g bands are represented by a TB Hamiltonian with diagonal elements of the 
form 

Hxy,xy{k) = ed + to(Cx + Cy) + hCxCy + [t2 + ts^Cx + Cy) + UCxCy\Cz, (l) 

where q = 2cos(A;ja) [i = x, y, z) and a is the lattice constant fll^. Cyclic permutations 
yield Hyz^yz{k) and Hzx,zx{k). The ti denote effective hopping integrals arising from the V- 
0-V hybridization. to,2, ^i,3, and specify the interaction between first, second, and third 
neighbors, respectively, and Itol ^ l^il, 1^21, l^sl, |^4|- Off-diagonal elements which vanish at 
high-symmetry points are small and are neglected here. 

In the case of LSFO and LSMO, we performed TB band-structure calculations explicitly 
including oxygen 2p orbitals for the three-dimensional perovskite structure in the antiferro- 
magnetic (AF) and ferromagnetic (FM) states. TB calculations for the paramagnetic (PM) 
state were already reported in Refs. Iisl . 14 1. The size of the matrix for the latter cal- 



culation was 14 X 14 (TM 3d orbitals: 5, oxygen 2p orbitals: 3x3 = 9). As for LSFO, 
we performed the calculation assuming the G-type AF state, where all nearest-neighbor Fe 
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atoms have antiparallel spins. The effect of the G-type antiferromagnetism was treated phe- 
nomenologically by assuming an energy difference AE between the spin-up and spin-down 
Fe sites. As for LSMO, we performed the calculation assuming the FM state. The effect 
of ferromagnetism was again treated phenomenologically by assuming an energy difference 
AE between the spin-up and spin-down Mn sites. The size of the matrix to be diagonal- 
ized was 28 x 28 for LSFO (G-type AF), and 14 x 14 for LSMO (FM). Parameters to be 
adjusted are AE, the energy difference between the TM 3d level and the O 2p level — e^, 
and Slater-Koster parameters (pda), (pdn), (ppa), and (ppvr) Here, is the average 

of the majority-spin d-\eve\ — AE/2 and the minority-spin d-level -|- AE/2. The ratio 
{pda)/ (pdn) was fixed at —2.2 and (ppcr) and (ppTf) were fixed at 0.60 eV and —0.15 

eV, respectively, (pda) is expected to be in the range of —(1.4 — 1.9) eV for LSFO and near 
—1.8 eV for LSMO from a configuration-interaction (CI) cluster-model calculation reported 
in Refs. |3, 3- Crystal-field splitting lODq of 0.41 eV was taken from Ref. 1^. Crystal- 
field splitting of O 2p orbitals is not taken into account. The TB basis functions used in the 
calculations of LSFO and LSMO are given in Tables HI and IIIH and the matrix elements of 
the Hamiltonians for LSFO and LSMO are given in Tables HTl and HV| respectively. 

Figure H shows the Brillouin zones (BZs) for the perovskite structure in the PM and FM 
states [(a)], the G-type AF state [(b)], and the G-type AF state with GdFeOs-type distortion 
[(c)]. The BZ of the FM state is the same as that of the PM state. The BZ of the G-type 
AF state is that of the face-centered cubic (fee) lattice. In this paper, we will neglect the 
effect of the GdFeOs-type distortion and will only consider the effect of magnetic ordering. 

Figure 121 shows schematic band structures along the F-X line obtained by TB calculation. 
Here, we have fixed — Cp — AE/2 at 3 eV. Panel (a) shows the band structure of the PM 
state calculated including only p-p interactions, namely, without p-d interaction {{pda) = 0). 
There are significant dispersions in the O 2p bands caused by the p-p interaction, whereas 
there are no dispersions in the Cg and t2g bands, where there exists only the crystal-field 
splitting of lODq. Panel (b) shows the band structure of the PM state calculated with 
{pda) = —2.0 eV. The eg and t2g bands now become dispersive due to the p-d interaction. 
Panels (c) and (d) show the band structure of the FM and G-type AF states, respectively. 
The values of AE were taken as 5.0 eV in both cases. In the FM state, since ed — ep — AE/2 is 
fixed, the up-spin (majority-spin) Cg and t2g bands do not change appreciably compared with 
those in the PM state, whereas the down-spin (minority-spin) Cg and t2gbands are pushed 
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FIG. 1: Brillouin zones of the perovskite structure. (a)PM and FM states. (b)G-type AF state. 
Broken lines represent the BZ of the PM state. (c)G-type AF state with GdFeOa-type distor- 
tion. Broken hnes and dotted hues represent the BZs of the PM state and the G-type AF state, 
respectively. 



up by AE (out of the range of Fig. 12 (c)). In the G-type AF state, band dispersions along 
the R-M line are folded onto the F-X line, and hybridization between these bands causes a 
maximum in the Cg bands in the middle of F and X. 



kz broadening 

In ARPES measurements, there is an intrinsic uncertainty in the momentum of electrons 
in the solid perpendicular to the crystal surface due to a finite escape depth A of photoelec- 
trons Q|. Here, the 2;-axis is taken as the direction perpendicular to the surface toward the 
vacuum side. Inside the solid {z < 0) the wave function of an emitted photoelectron with a 
wave number kzo is approximated by 

ip{z) = -^exp{ikzoz)exp{z/2X), (2) 

where A denotes the escape depth of a photoelectron. = exp(2;/A)/A shows an 

exponential decay of photoelectrons within the solid. The Fourier transform of Eq. Q, 
(pih), is 

and 10(^2) p becomes proportional to 

^^^^^ = 2^{kz-kzor+ii/2xr 
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FIG. 2: Band structures obtained by TB calculation, — — A£'/2 = 3 eV has been fixed. 
(a)PM state with only p-p interaction included {{pda) = 0). (b)PM state with {pda) = —2.0 eV. 
(c)FM state with AE = 5.0 eV. (d)G-type AF state with AE = 5.0 eV. 

which means that the intrinsic uncertainty of kz is given by the Lorentzian function with 
the full width at half maximum of 1/A. 

The effect of the kz broadening is very important for the interpretation of ARPES results. 
For example, we consider an isotropic band with the bottom at the T point, expressed as 

E^a{kl + kl + kl). (5) 

If we measure photoelectrons with momentum 1^ — {K^-Ky, K^), where Kz is related to kzo 
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through the momentum conservation kx = and ky = Ky and the energy conservation, 
the obtained spectrum becomes, 



/oo 
dk,L{kz)S{E - a{kx^ + ky" + A;/)), 
-oo 

From Eq. Q and using the identity 



(6) 



5{x — akz^) 



2\/ax 



5[k,-J-\ + S(k, + J- 



a 



we obtain 



D{kxi ky, E) 



1 



a 



:X 



(7) 



47rX^a{E ~a{k/ + ky')} 
1 



i^^{E-a{kJ + ky')}/a-k,o) +{l/2Xy 



(8) 



+ 



{E-a{kJ + ky')}/a + k,o) +(1/2A)2 



We plot the function D{kx, ky, E) with k^ = ky = for various values of A in Fig. 01 Here, we 
have assumed that a = 8 eV-A , which corresponds to m* ~ O.brrie {m* is an effective band 
mass and me is the electron mass in vacuum), and fc^o = 7r/4c, where c = 3.830 A, the out- 
of-plane lattice constant of LSMO (x = 0.4) thin films. One can see that in bulk- sensitive 
measurements using high-energy photons, like A = 20 A, there is a peak at E = ak^o' . In 
contrast, in surface- sensitive measurements using low-energy photons, there is no peak at 
akzo' but D(E) is peaked around ~ 0. The density of states at the F point is very large 
because the band has a van Hove singularity at this point. Since bands usually have van 
Hove singularities along a high-symmetry line, the effect of k^ broadening is very important 
when we trace near the high-symmetry line using low energy photons. Our photon energies 
used in the ARPES measurements were in the range of 20 - 100 eV and therefore A ~ SA, 
meaning that considering the fc^-broadening effect is crucial in the interpretation of ARPES 
results. 



7 





• 

■' 

■ ; 

\ / 


1 ■\ 
\ X = 1 A 

> OA 

\ X = 5 A 

■• X = 1 A 

X = 20 A 

\ \ 

V. 

\ 

■ — ■ ^ > 




1 =Ei 



akj 0.5 1.0 



Relative energy (eV) 

FIG. 3: ARPES spectrum D{kx,ky, E) with kx = ky = for various values of escape depths of 
photoelectrons A. 



RESULTS AND DISCUSSION 



SrVO. 



SVO is an ideal system to study the fundamental physics of electron correlation because 
it is a prototypical Mott-Hubbard-type with electron configuration. For the bandwidth- 
control system Ca1__i.Sr2.VO3 (CSVO), the existence and the extent of effective mass en- 
hancement and spectral weight transfer have been quite controversial. Early photoemission 
results have shown that, with decreasing x, i.e., with decreasing bandwidth, spectral weight 
is transferred from the coherent part to the incoherent part [2^ in a more dramatic way than 
the enhancement of the electronic specific heats 7 |2^. In contrast, in a recent bulk-sensitive 
photoemission study using soft x rays it was claimed that there was no appreciable spectral 



weight transfer between SrVOa and CaVOg 



221. 



ARPES measurements on SVO were performed at beamline 5-4 of Stanford Synchrotron 
Radiation Laboratory (SSRL). Bulk single crystals of SVO were grown using the traveling- 
solvant floating zone method. Samples were first aligned ex situ using Laue diffraction, 
cleaved in situ along the cubic (100) surface at a temperature of 15 K and measured at the 
same temperature. Details of the experimental conditions of ARPES are described in Ref. 

An example of ARPES spectra for SVO near Ep are shown in Fig. |3J The coherent part 
within ~ 0.7 eV of Ep shows a parabolic energy dispersion, consistent with the view that 
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FIG. 4: ARPES spectra of SrVOs along momentum cut 2 in the Brillouin zone shown in Fig. 
|5J (a) Intensity plot in the E-ky plane. Dispersive feature within ~0.7 eV of Ep is the coherent 
part, while a broad feature around -1.5 eV is the incoherent part, (b) Energy distribution curves 
(EDC's). 
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FIG. 5: Band structure and the Fermi surface of SrVOs. (a) Band dispersions obtained by local- 
density approximation (LDA) band-structure calculation 23]. (h)ky=0 cross-sectional view of the 
Fermi surfaces and the momentum loci corresponding to the mapping in Fig. 



the coherent part corresponds to the calculated band structure. On the other hand, the 
incoherent part centered at ~ -1.5 eV does not show an appreciable dispersion but shows a 
slight modulation of the intensity, which coincides with that of the coherent part, as seen in 
FiglU^a). To compare with the band-structure calculation, we shall focus on the dispersive 
feature of the coherent part. 
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FIG. 6: Spectral weight mapping at Ep for SrVOs. (a) Mapping with hiy=28.5eV. (b) Mapping 
with hiy= 24 eV. Note that the mapping is projection on the kx-ky plane. 

The band dispersions obtained by local-density approximation (LDA) band-calculation 
23[ are shown in Fig. (a). The hopping parameters of the TB model (1) have been 
obtained by fitting the TB bands to the results of LDA band-structure calculation [2^. 
Figure El (b) illustrates the ky=0 cross-sectional view of the Fermi surfaces expressed by the 
TB model. The three bands correspond to the dxy, dyz and dzx orbitals of V and scarcely 
hybridize with each other. Therefore, these dispersions are almost two-dimensional and 
form three cylindrical Fermi surfaces intersecting each other around the F point, although 
the crystal structure is three-dimensional. The loci of electron momenta in the kx-k^ plane 
for constant photon energies hv= 24 eV and 28.5 eV are plotted in Fig. El^b). Here, the 
inner potential of 10 eV has been assumed. 

Figure IHl (a) and (b) shows spectral weight mapping in space at Ep on the 

momentum loci for the photon energies /iz/=28.5 and 24 eV, respectively. The spectral 
weight for both photon energies are largely confined within the cylindrical Fermi surfaces 
extended in the kz direction. According to Figl^l^b), cut 2 corresponds to momenta near 
the X point [k = (1,0,0)]. Therefore, the observed spectral weight comes only from one 
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FIG. 7: ARPES spectra of SrVOa for cuts 1 and 2 in Fig. [51 Panels (a) and (b) show intensity 
plots in E-k space for cuts 1 and cut 2, respectively. Angle-independent backgrounds have been 
subtracted. Calculated TB band are also superimposed. Circles in (b) are peak positions of MDC's 
and represents band dispersion. Panels (c) and (d) are EDC's corresponding to panels (a) and (b), 
respectively. Vertical bars are guides to the eye indicating the positions of the dispersive features. 



cylindrical Fermi surface referred to as the 7 sheet j2J], which arises from the d^y orbital. 
On the other hand, the spectral weight in the first BZ for hu =28.5 eV (cut 1) is enclosed 
by the three cylindrical Fermi surfaces [2^. Indeed, EDC's along cut 1 [Fig. [3 (c)] show 
complicated features suggestive of two energy dispersions. 

Figure Wi^) and (b) show intensity plots in the energy-momentum [E-k) space corre- 
sponding to each cut in Fig. |S| The E-ky plots are compared with the energy dispersions 
expected from the LDA calculation Q]. Figure [7||^c) and (d) show EDC's corresponding to 
panels (a) and (b), respectively. The EDC's of cut 1 in panel (c) show parabolic energy 
dispersions with the bottom at ~-0.5 eV, shallower than that expected from the LDA cal- 
culation ~ -0.9 eV. Since the momenta of cut 1 are within the a and (3 sheets of the Fermi 
surface, there are two dispersions. These dispersions may correspond to the dispersions of 
the dxy and dy^ bands. The dispersion of d^x is not clearly observed in the present plot, 
probably due to the effect of transition matrix-elements. 

For cut 2 [Fig. [7| (b) and (d)], since it does not intersect the dyz and dzx bands, only 
the dxy band are observed. In panel (b), we have derived the band dispersion for cut 2 
from MDC peak positions indicated by circles. By comparing the Fermi velocity of the LDA 
band structure and that of the present experiment, one can see an overall band narrowing 
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in the measured band dispersion. Near Ep, we obtain the mass enhancement factor of 
m* /nib ~ 1-8 ± 0.2, which is close to the value m*/mb = 1.98 obtained from the specific 
heat coefficient 7 j21|. Since the LDA calculation predicted that the d^y band is highly 
two-dimensional and nearly isotropic within the k^-ky plane, a similar mass enhancement is 
expected on the entire Fermi surface. 

It has been pointed out that the photoemission spectra of SVO taken at low photon 
energies have a significant amount of contributions from surface states, particularly in the 



incoherent part 



22 



25|. In order to distinguish bulk from surface contributions in the 



present spectra, we consider the character of surface states created by the discontinuity of 
the potential at the surfaces. According to a tight-binding Green's function formalism for a 
semi-infinite system 11, 2^, surface-projected density of states of surface-parallel d^y band 
and surface-perpendicular dy^^z^ band have different nature. The surface dxy band does not 
show an appreciable change from the bulk band because it has the two-dimensional character 
within the pl«>e. In contrast, =pect.al w«ght of the surface states . .edistrfbuted 

between the bulk dyz,zx band and the Fermi level p|| , because kz is no longer a good quantum 
number. This leads to the absence of clear dispersions of the surface dyz^zx bands in ARPES 
as demonstrated in Ref. 0. Accordingly, the observed dispersion of the d^y band should 
represent that of the bulk d^y band. 



Lai_a;Sra;Fe03 

LSFO undergoes a pronounced charge disproportionation and an associated MIT around 
X ~ 2/3 [27]. One striking feature of LSFO is that the insulating phase is unusually wide 
in the phase diagram (0 < a; < 0.5 at room temperature and even < a; < 0.7 at low 
temperatures) 28|. In a previous photoemission study, the gap at the Fermi level (Ep) was 
observed for all compositions for < x < 0.67 19j, consistent with the wide insulating 
region of this system. 

ARPES measurements on LSFO were carried out using a photoemission spectroscopy 
(PES) system combined with a laser molecular beam epitaxy (MBE) chamber, which was 
installed at beamline BL-IC of Photon Factory, KEK The preparation and character- 
ization of the LSFO films are described in Ref. By low energy electron diffraction 
(LEED), sharp 1x1 spots were observed with no sign of surface reconstruction, confirming 
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the well-ordered surfaces of the films. Details of the experimental conditions are described 
in Ref. Q. 

The left panel of Fig. |S1 (a) shows a gray-scale plot of the experimental band structure 
obtained by ARPES measurements of an LSFO {x = 0.4) thin film with the photon energy of 
74 eV Isi]. Here, the second derivatives of the EDC's are plotted on the gray scale, where dark 
parts correspond to energy bands. In Fig. |H1 (b), hole pockets obtained by a TB calculation 
(described below) are also shown. The trace obtained for hu = 74 eV crosses the calculated 
hole pocket. The structure at —1.3 eV shows a significant dispersion, and is assigned to 
the Fe 3d majority-spin eg bands. The structure at —2.4 eV shows a weaker dispersion 
than the Cg bands, and is assigned to the Fe 3d t2g majority-spin bands. The structures at 
— (4 — 7) eV are assigned to the O 2p bands. The Cg bands do not cross Ep, consistent with 
the insulating behavior and the persistence of the gap observed in the angle-integrated PES 
(AIPES) spectra 



' Fjg. IHl (a)] to the observed band dispersions 
has been obtained with e^- = 2.65 eV, 



The best fit of TB calculation [right panel o: 
and the optical gap of 2.1 eV of LaFeOs (LFO) 
{pda) = —1.5 eV, and AE = 5.3 eV. We also considered the effect of kz broadening given 
by a Lorentzian function as described in Sec. 2.2. Here we have used the value A = 5 A. 
Akz ~ 1/A (~ 0.2 A ^) is approximately 10 % of the Brillouin zone (27r/a). Here, it should 
be noted that — Cp, (pda) and AE primarily determine the Fe 3d — O 2p band positions, 
their dispersions, and the optical band gap, respectively. The dispersions of the Cg bands 
were successfully reproduced. The very weak dispersions of the t2g bands and the width 
of the O 2p bands were also well reproduced by this calculation. The energy position of 
the calculated Ep determined from the filling of electrons was not in agreement with the 
experimental Ep, however. This discrepancy corresponds to the fact that this material is 
insulating up to 70 % hole doping while the rigid-band model based on the band-structure 
calculation gives the metallic state. 

Lai :,.Sr:,.Mn03 

LSMO has been attracting great interest because of its intriguing properties, such as 
the CMR effect Q 

is an antiferromagnetic insulator, and hole-doping induced by the substitution of Sr for La 
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FIG. 8: Comparison of ARPES spectra of an LSFO thin film taken at hu = 74 eV with TB 
calculation. The left panel of (a) shows the experimental band structure (dark parts correspond to 
energy bands) and the right panel of (a) shows the result of TB calculation. The traces in /c-space 
and calculated hole pockets are shown in (b). 



produces a ferromagnetic metallic phase (0.17 < x < 0.5) 

ARPES measurements on LSMO were carried out using the same system as in the case of 
LSFO. The preparation and characterization of the LSMO films are described in Refs. js^. 
The surface structure and cleanliness of the films were checked by LEED. Some surface 
reconstruction-derived spots were observed in addition to sharp 1x1 spots, confirming 
the well-ordered surfaces of the films, but surface reconstruction is expected to affect the 

n 

obtained band dispersions. Details of the experimental conditions are described in Ref. j8||. 
The left panel of Fig. El (a) shows a gray-scale plot of the experimental band structure 



obtained 
of 88 eV 



jy ARPES measurements of an LSMO (x = 0.4) thin film with the photon energy 
. Here, the second derivatives of the EDO's are plotted on the gray scale, where 



dark parts correspond to energy bands as in the case of LSFO along the trace in fc-space 
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FIG. 9: Comparison of ARPES spectra of an LSMO thin film taken at hv = 88 eV and TB 
calculation. The left panel of (a) shows the experimental band structure (dark parts correspond 
to energy bands) 8] and the right panel of (a) shows the result of TB calculation. The traces in 
A;-space and calculated electron pockets are shown in (b). 



shown in Fig. El (b). Note that this trace is approximately along the F-X direction. We shall 
assume that the trace is exactly along the F-X direction for the moment. The effect of the 
deviation from the F-X line will be discussed below. Here, electron pockets obtained by a 
TB calculation are also shown in Fig. IHl (b). One of the Mn 3d majority-spin eg bands at 
— (0 — 1) eV shows a significant dispersion and crosses Ep, corresponding to the half metallic 
behavior The Ep crossing of the Cg bands is consistent with the observed Fermi edge 
in the AIPES spectra js^. The Mn 3d majority-spin bands at —(1 — 2) eV show weaker 
dispersions than the eg bands. The bands at —(2 — 3) eV shows dispersions with a shorter 
period in fc-space than the Brillouin zone. This means that these bands are affected by 
the surface reconstruction detected in the LEED pattern. The O 2p bands are located at 
-(4-7) eV. 

The best fit of the TB calculation to the observed band dispersions [right panel of Fig. 1^1 
(a)] has been obtained with — = 5.3 eV, (pda) = —2.0 eV, and AE = 4.6 eV. The 
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strong dispersions of the Cg bands, the weak dispersions of the t2g bands, and the width of 
the O 2p bands were successfully reproduced, although there are no bands in the calculation 
corresponding to the experimental surface bands. The bottom of the Cg bands is at ~ —0.5 
eV in the experiment and at ~ —1.5 eV in the calculation. This difference may come from 
the mass renormalization caused by strong electron correlation in this system. One can 
estimate the mass enhencement to be m*/mb ~ 2.5. 

Figure Uni (a) shows the calculated density of states (DOS) of the FM state. The partial 
DOS's for the majority-spin Mn Cg, the majority-spin Mn t2g, the minority-spin Mn Cg, and 
the minority-spin Mn t2g orbitals are shown in the lower panels. No band gap opens between 
the majority-spin Cg bands and the minority-spin t2g bands for the present AE value. In 
Fig. ^1 (a), the Ep position has been determined from the filling of electrons, and also 
calculated DOS at Ep [D^{Ep)] has been obtained. Using the relationship 

7" = ^kB'D\Ep), (9) 

we obtain the calculated electronic specific heat coefficient by 7^ ~ 1 [mJ mol^ K~^]. The 
experimental electronic specific heat coefficient has been reported as 7 ~ 3.5 [mJ mol^ K~^] 
33[. Thus m* /rrih = 7/7?, ~ 3.5, in fairly good agreement with m* /rrih ~ 2.5 derived from 
ARPES, which means that effective mass m* is renormalized as 

m*/mb ~ 2.5 - 3.5. (10) 



Figure ITUl fb) shows the comparison of the chemical 



potential shift obtained from the TB 



30 



within the rigid-band model 



calculation in the FM region of LSMO (0.2 < x < 0.4 
(i.e., assuming the constant exch ang e splitting) and that determined experimentally from 
core-level photoemission spectra j3^. The experimental chemical potential shift is slower 



than the calculated one by a factor of ~ 0.4, which means 

dfi 0.4 

dn ^ D\Epy 

For a metallic system, the rigid-band picture predicts that 

d^i l + F fmh\ l + F 



(11) 



(12) 



dn D*{Ep) \m*J D^EfY 

where D*{Ep) is the DOS of renormalized quasiparticles (QPs) and F is a parameter which 
represents the effective QP-QP repulsion. From Eqs. (fTU Il -lfT ^ . one obtains F ~ — 0.4. 

16 




" -0.2 



-10 -5 5 

Energy (eV) 



10 



0.0 0.2 0.4 0.6 
Sr Concentration x 



FIG. 10: Chemical potential shift in LSMO. (a)DOS obtained from TB calculation. (b)Comparison 
of the chemical potential shift calculated for the rigid-band-model using TB calculation and that 
determined experimentally from core-level photoemission spectra js^. 



This value is much smaller than the values of Fq in AF compounds. For example, it was 
reported that Fq ~ 7 for La2_-a;Sra;Cu04, and ~ 6 for Lai-^Sr^VOs 3^. In the FM case, 
one obtains small values of F compared to the AF case, probably because of the strong 
exchange interaction between carriers in the FM case. 

Now we go back to the effect of the deviation of the trace with hu = 88 eV from the F-X 
direction. Figure ITD shows a TB simulation of the ARPES spectra of LSMO taken at 88 eV 
using the above-mentioned parameter values. Panels (a), (b) and (c) show the results for 
Akz = (without kz broadening), Ak^ = 1/A, and Ak^ = oo (uniform integration in k^), 
respectively. In panel (a), the energy position of the bottom of the Cg bands at = A 
is different from that at A;|| = 1.6 A ^, while in panel (b) the bottoms are at almost the same 
energy positions. In experiment [Fig. IHl (a)], the bottoms are almost at the same energy 
positions, which means that the experimental results are in good agreement with (b), that 
is Akz = 1/A. In panel (c), there are two bands approaching Ep, in disagreement with 
experiment. 

Figure ^1 shows comparison of the dispersion along the F-X direction (dots) and simu- 
lated spectra (gray scale) obtained assuming Akz = 1/A for hu = 88 eV. There are some 
discrepancies in the high-binding-energy region —(4 — 7) eV, but the overall agreement is 
rather good in the region near Ep [—(0 — 2) eV]. This means that one can trace the band 
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FIG. 11: Simulated ARPES spectra of LSMO {x = 0.4) taken at 88 eV on TB calculation, (a) 

Ak^ = 0, (b) Ak^ = 1/A, (c) Ak^ = oo. 



dispersions almost along the F-X direction by performing ARPES with a fixed photon en- 
ergy of hu = 88 eV because of the large fc^-integrated DOS at a high-symmetry line (in 
this case the F-X direction). Even if the trace somewhat deviates from a high-symmetry 
line, one can trace the dispersion along the high-symmetry line due to the large DOS and 
the effect of kz broadening. By changing the photon energy with changing the emission 
angle in ARPES measurements, one can map band dispersions precisely along the F-X line. 
However, the photoionization cross section varies with a photon energy, making the interpre- 
tation difficult. In practice, fixing a photon energy and changing an emission angle is a more 
efficient and useful way to obtain band dispersions in such materials. We also revealed that 
considering the effect of kz broadening is important for the interpretation of the ARPES of 
the three-dimensional materials. 



CONCLUSION 

We have analyzed the ARPES data of SVO bulk crystals and LSFO {x = 0.4) and LSMO 
(x = 0.4) thin films using TB band-structure calculations. As for SVO, the enhanced 

effective electron mass obtained from the energy band near Ep was consistent with the 
bulk thermodynamic properties. The experimental band structures of LSFO and LSMO 
were interpreted using a TB band-structure calculation. The overall agreement between 
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FIG. 12: Comparison of the dispersion along the F-X direction (dots) and simulated spectra (gray 
scale) for LSMO {x = 0.4) calculated assuming Akz = 1/A at 88 eV. 

experiment and calculation is fairly good. However, in the case of LSFO, the energy position 
of the calculated Ep was not in agreement with the experimental Ep, reflecting the insulating 
behavior. In the case of LSMO, experimental energy bands were narrowed compared with 
calculation due to strong electron correlation. We also demonstrated that considering the 
effect of kz broadening is important for the interpretation of the ARPES results of such 
three-dimensional materials. 
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BASIS FUNCTIONS AND MATRIX ELEMENTS FOR TB BAND-STRUCTURE 
CALCULATIONS 



TABLE L TB basis functions used in the band-structure cal- 
culation of LSFO 





No. 


Origin 


Function No. 


Origin 


Function 


Fe 3d 


1 


(0, 0, 0) 


xy 


15 


0(1,0, 0) 


xy 




2 


(0,0,0) 


yz 


16 


0(1,0,0) 


yz 




3 


(0,0,0) 


zx 


17 


a(l,0,0) 


zx 




4 


(0,0,0) 


3z^ - 


18 


0(1,0,0) 


3z^ - 




5 


(0,0,0) 


x^ — y^ 


19 


0(1,0,0) 


9 9 

— 


2p 


6 


50(1, 0,0) 


X 


20 


^a(3,0,0) 


X 




7 


ia(l,0,0) 


y 


21 


ia(3,0,0) 


y 




8 


ia(l,0,0) 


z 


22 


ia(3,0,0) 


z 




9 


ia(0,l,0) 


X 


23 


ia(2,l,0) 


X 




10 


50(0,1,0) 


y 


24 


ia(2,l,0) 


y 




11 


50(0,1,0) 


z 


25 


50(2,1,0) 


z 




12 


ia(0,0, 1) 


X 


26 


ia(2,0,l) 


X 




13 


50(0,0,1) 


y 


27 


ia(2,0,l) 


y 




14 


50(0,0,1) 


z 


28 


ia(2,0,l) 


z 



TABLE II: Nonzero matrix elements of the TB Hamiltonian 
used in the calculation of LSFO. ^ = k^a, 77 = kya, and 
C = kza. 

A. d-d interactions 

Hi^i = H2,2 = Hs^s = ed- ^Dq - AE/2, i/4,4 = i^5,5 = erf + QDq - AE/2 
Hi5,i5 = i?i6,i6 = ^^17,17 = ed- 4.Dq + AE/2, His,i8 = i?i9,i9 = + ^Dq + AE/2 
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B. p-p interactions 

Ci = \{pp(y) + \{vpT^), C2 = \{vp<y) - \{vpT^), C3 = ivpT^) 

H^fi = i?10,10 = -H^14,14 = -f^20,20 = -^24,24 = -^28,28 = 
H7J = i?8,8 = ^9,9 = ^^^11, 11 = -ffl2,12 = -ff 13,13 = 
-f^21,21 = -f^22,22 = ^^23,23 = -^25,25 = -^26,26 = -^27,27 = 
-f^6,9 = -^7,10 = -^20,23 = -^21,24 = 2Ci COS (^^^ - ^-qj 
He,10 = -^7,9 = -^20,24 = -f^21,23 = — 26*2 COS (^^^ - ^7]^ 

He,i2 = Hs^u = i!f20,26 = -^22,28 = 2Ci cos (^^^ — ^c) 

-f^6,14 = -f^8,12 = -^20,28 = -f^22,26 = — 26*2 COS — 
^^6,23 = ^^7,24 = -^20,9 = -f^21,10 = 2Ci COS + i??) 

-f^6,24 = -f^7,23 = -^20,10 = -f^21,9 = 2C2 COS (^^ + ^7]^ 
Hq,26 = -^8,28 = -^20,12 = -f^22,14 = 2Ci COS (^^^ + 
-f^6,28 = -^8,26 = -f^20,14 = -^22,12 = 26*2 COS (^^^ + 
i?7,13 = ^21,27 = 2C3 COS (i^ - 
i?7,27 = ^21,13 = 2C3 COS + ^C) 
^^8,11 = i?22,25 = 2C3 COS (i^ - ir/) 
^8,25 =^^22,11 = 2C3 COS (i^ + ir/) 
-^9,12 = -ff23,26 = 2C3 COS (^^T] - ^(j 
H9,26 = i?23,12 = 2C3 COS (^7/ + 
-f^l0,13 = -ffll,14 = -^24,27 = -^25,28 = 2Ci COS (^^Tj — 
-^10,14 = -f^ll,13 = H24,28 = -f^25,27 = — 2C2 COS (^^7] — 
-^10,27 = -ffll,28 = -f^24,13 = -f^25,14 = 2Ci COS (^^7] + i^) 
-f^l0,28 = -ffll,27 = -f^^24,14 = -f^^25,13 = 2C2 COS (^^7? + ^(j 



C. p-d interactions 

Hi,7 = i?3,8 = -H"i5,2i = -ffi7,22 = -{pdir) exp (il^^ 
Hi,9 = H2,u = -H'15,23 = -f^i6,25 = -{pdir) exp {i^'n'j 
Hi,2i = H^,22 = Hi5j = Hn^s = (pdir) exp (-i^C) 

-H"l,23 = -f^^2,25 = ^^15,9 = ^^16,11 = {pdir) exp (-^5??] 

^^2,13 = -^3,12 = -H"i6,27 = -H"i7,26 = -{pdn) exp 
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TT 

— -"3,26 


TT 

— -"16,13 - 


= -ffi7,i2 = {pdir) exp 1 


TT 

-"4,6 


TT 

— -"18,20 


= l{pda)ex.p (il^^ 


jj 

-"4,10 


u 

— -"18,24 


= ^[pda) exp (^i^Tj^ 


TT 


TT 

= -"18,28 


= -{pda)exp (i^C^ 


TT 

-"4,20 = 


TT 

= -"18,6 = 


-2(pdcr)exp [-12^ 


TT 

-"4,24 = 


TT 

- -"18,10 = 


— ^{pda) exp (^—i^V^ 


TT 

-"4,28 


TT 

— -"18,14 


= {pda)exp 


^^5,6 = 


-f^l9,20 = 


-^(pda)exp(iic) 


^^5,10 ■ 


= -H^19,24 = 


= ^(pcZcr)exp (^^5?7^ 


^^5,20 = 


= ^19,6 = 


^(pd(7)exp(-zic) 


-f^^5,24 = 


-f^l9,10 = 


— ^{pda) exp (— i^r/ 



TABLE III: TB basis functions used in the band-structure 
calculation of LSMO 



No. Origin Function No. Origin Function 



1 


(0,0,0) xy 


2p 6 


ia(l,0,0) 


X 


2 


(0,0,0) yz 


7 


ia(l,0,0) 


y 


3 


(0, 0, 0) zx 


8 


ia(l,0,0) 


z 


4 


(0,0,0) - 


9 


ia(0,l,0) 


X 


5 


(0,0,0) x2-y2 


10 


ia(0,l,0) 


y 






11 


ia(0,l,0) 


z 






12 


ia(0,0,l) 


X 






13 


|a(0,0,l) 


y 






14 


ia(0,0,l) 


z 
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TABLE IV: Nonzero matrix elements of the TB Hamiltonian 

used in the calculation of LSMO. ^ = kxa, rj = kyU, and 
C = kza. 

A. d-d interactions 

= H2,2 = -^3,3 = ed- 4:Dq - AE/2 (up spin), - 4Dq + AE/2 (down spin) 
H4,4 = -ff5,5 = ed + QDq - AE/2 (up spin), + QDq + AE/2 (down spin) 

B. p-p interactions 

Ci = hiPP^) + Upp^)^ C's = Upp^) - Upp^)^ C's = (pptt) 

^6,6 = -^10,10 = -ffl4,14 = 

H7J = Hs^s = -f^9,9 = = -f^l2,12 = -H^13,13 = 

^^6,9 = Hr^w = 4Ci cos(iO cos(ir/) 
He,io = H7^9 = -4^2 sin(ie) sm{^ri) 
H%,i2 = H8,u = 4Ci cos(^0 cos(^C) 
^^6,14 = Hs,i2 = -4C2 sin(ie) sin(iC) 
H7,i3 = cos(iO cos(iC), ^^8,11 = 4C3 cos(iO cos(ir?) 
-^9,12 = 4C3 cos(ir/) cos(iC) 
-H"io,i3 = -H"ii,i4 = 4Ci cos(i?7) cos(iC) 
i^io,i4 = -^"11,13 - 4C2 sin(577) sin(^C) 

C. p-d interactions 

Hi,7 = H3,s = -2i{pd7r) sin(iO 
Hi,9 = H2,u = -2i{pd'K) sm{\r]) 
H2,i3 = H^,i2 = -2i{pdTr) sin(iC) 
H4,6 = i{pda) sin(i^), = i{pda) sm{^ri) 

if4,i4 = -2i{pda) s'm{^C), ^^5,6 = -\/3i(pd(j) sin(i^) 
H5,io = V3i{pda) sin(ir/) 
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